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Abstract 

This report analyzes the eddy currents induced in a solid conducting 
sphere by a sinusoidal current in a circular loop. Analytical expressions for 
the eddy currents are derived as a power series in the vectorial displacement 
of the center of the sphere from the axis of the loop. These are used for 
first order calculations of the power dissipated in the sphere and the force 
and torque exerted on the sphere by the electromagnetic field of the loop. 
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Table of Notation 


E -> Electric Field 

B -* Magnetic Induction 

A -» Vector Potential 

J -> Current Density 

I -» Current 

<•> -» Frequency 

* 

// o ■+ Permeability of free space 

5 -♦ Skin depth 

a -* Conductivity 

-+ Unit Radial Vector 

-* unit Latitudinal Vector 

■+ Unit Azimuthal Vector 

e x ,e y ,e x -+ Unit Cartesian Vectors 

s-(s,6 s ,* s ), o(c,© c ,* c ), r-(r,0,*) -» Postion vectors 

d -> Displacement of sphere center from axis of loop 
a -» Radius of sphere 

b -» Radius of loop 

n -* Normal unit vector to loop 
P n ,P^ -> Legendre functions 
J ,K ,1 ,Y -* Bessel functions 
C ,A ,C ,A ,B ,D -* constants 
-» Kronecker delta 

i -» Unit imaginary 

Xn m •* Vector Spherical Harmonics 
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Introduction 


This report is an extension of research by David Sonnabend [1,2] on the 

control of suspended objects by eddy current forces, and it was carried out at 

his suggestion. The results are potentially useful in the design of a quasi 

% 

drag free gradiometer to be used in artificial space satellites. 

In this study we are concerned with effects on a conducting sphere of the 
electrooagnetic field produced by a steady state alternating current in a 
circular coil. The oscillation field induces eddy currents in the conductor 
which dissipate power and interact with the field to produce a force and 
torque on the sphere. Our objective is to derive analytical expressions 
describing these effects. The central task is to solve the steady state 
electromagnetic boundary value problem for a conducting sphere in the 
oscillating field of the current loop. The problem was solved some time ago 
for the special case where the center of the sphere lies on the axis of the 
loop [3,4]. Tegopoulos and Kriezis [4] give a valuable survey of analytical 
work on problems of this type. Hannakam [5] has found an elegant closed 
solution to our general problem in terms of an integral over the loop. His 
result is summarized in Appendix M of this report. However, after studying 
the difficulties in evaj.ucu..i.iiy Ills 4 ‘ Ka4 ‘ urn,1 d be more 

practical to solve the problem ab initio by a different method. We studied a 
variety of promising approaches. Expansions in terms of Lame' polynomials [6] 
seemed most promising because these functions possess both the symmetries and 
asymmetries of our problem. Unfortunately, the mathematical theory of Lame' 
polynomials is not sufficiently well developed to make all the necessary 
calculations easy. We finally settled on expansions in terms of spherical 
harmonics in large measure because many theorems about these functions are 
available to facilitate calculations. This enabled us to find a practical (if 
not optimal) solution to our problem. 
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To make our boundary value problem well defined and analytically 
tractable, we adopted the following idealizations of the physical situation: 

1) the coil is replaced by a single, circular current loop. If the coil 
needs to be more accurately modelled, this can be accomplished by 
considering it as a series of current loops and supercomposing the 
individual solutions. 

2) The current in the loop is taken to vary sinusoidally in time. That is 

I - I cos cot - Re {I e iwt } 

O O 

I o ,to are real constants. 

3) The propagation of the fields is regarded as instantaneous, that is, 
the fields change slowly with respect to their propagation time to 
points within the domain of this problem. Thus, all points in the 
conductor "see" a field of the same phase at the same instant in time, 
and the displacement current can be neglected. 

4) The conductor is isotropic, homogeneous, and non-magnetic with a 
relative permeability of one. Hence, within it Ohm's law (J * oE) 
applies and hysteresis can be ignored. 

5) The conductor is modelled as a solid sphere, ‘xms mooex vdiJ *. 

spherical shell as long as the skin depth is small relative to the 
thickness of the shell. 

6) The displacement of the center of the sphere from the axis of the loop 
is considered small compared to the radius of the sphere. This allows 
the use of an expansion about a position on the axis of the loop. 

7) Transcient effects are ignored and only steady state solutions are 
considered. 
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I. The electro-magnetic vector potential 

The first thing we do is introduce the idealizations of our problem into 
Maxwell's equations to derive a boundary value problem for the vector 
potential A. “Neglecting the displacement current, and considering ^materials 
with relative permeability of one, Maxwell's equations are: 

(1) V-E~p/e o 

(2) VxB - n o J 

dB 

(3) VxE - 

3t 

(4) V-B « 0 

From V*B - 0, it follows that B can be written as the curl of the vector 
potential A. 

(5) B - VxA 

The condition V*A - 0 can also be imposed to determine A uniquely. 
Substituting equation (5) into equation (3) yields 

d 3A 

(6) VxE (VxA) - V x ( ) 

3t 3A 

e<i wa ran falfp 

3A 

(7) E . 

3t 

Assuming dim’s law, J « oE and equation (7) yields 
3A 

(8) J - -a . 

3t 

We assume that the sources of the field are stationary and the field 
propagates instantaneously so it must have the same time dependence as the 
current loop. Accordingly, the vector potential can be written 
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A (r,t) - A (r) e iwt 

where only the real part is of physical interest. The eddy current density 
then is proportional to the vector potential 

(9) J - -iwcA. 

Substituting equations (5) and (9) into equation (2) we obtain 

(10) Vx(VxA) - i// o waA« k 2 A 

where k ■ ifJ o <aa. By virtue of the vector identity 
Vx( VxA) - V(V-A) - V 2 A 

and the relation V‘A-0, equation (10) reduces to Helmholz's equation 
(11a) V 2 A - k 2 A 

inside the sphere, and Laplace's equation 
(lib) V 2 A * 0 

immediately outside the sphere where the conductivity is zero. The solutions 
to equations (11a) and (lib) must be matched at the boundary of the sphere. 
According to Smythe [3, p. 305] the following boundary conditions must be 
satisfied. 

(1) A is continuous across the boundary 

(2) The normal component, of L 1- w W '.tir.u~_r 

(3) The tangential component of H ■ //B is continuous across any boundary. 
For this case v «= /j o inside the sphere as well as outside of it so the 
tangential componet of B is also continuous across the boundary. 

Thus, 

< 12 «> \l r .. - »J r .. 

(12b) VxAj... - t'xaj,... 

where A o is the total vector potential outside the sphere and A i is the total 
vector potential inside the sphere. 
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II. Calculating the vector potential of a current loop with respect to an 
origin off axis to the loop . 

To calculate the vector potential of a current loop, it is convenient to 
choose an origin on the symmetry axis of the loop. However, in calculating 
the total vector potential in the sphere, the symmetry of the sphere is most 
easily exploited by using an origin at its center. TO benefit from both 
choices of origin, the vector potential of the loop is first calculated with 
respect to an origin on the axis of the loop and then is expanded about the 
center of the sphere. This method would be inappropriate for large 
displacements of the center of the sphere from the axis of the loop. The 
geometrical parameters in our problem are shown in figure 1. 

The vector potential of a line current is found from 


(13) Aj 


4n J |R| 


/ 


where I is the current and can depend on time, 

dl « dl 1 is the line element directed along the current, 

|r| is the distance from dl to the field point. 

Extsandinq I Rp 1 in Legendre polynomials and integrating over the current 
loop we find (Appendix A), 


(14) A l 


// lb «° s n 

— i . 

2 n«l c n+ n(n+l) 


1 1 A 

P (n*c) P (n*s) *s 

n n 


|s|<|c| 

b is the radius of the loop 
s is the field point 
c is a vector to the loop 

p*(n*c), P*(n*s) are Associated Legendre functions 
*s is the azimuthal unit vector 
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Using s n P*(n*s)*s - V s x[P n (n*s)s n s], (Appendix B) 

H o Ib - P^n-S) A 

(15) A, - l — V x[p (n-s)s •]. 

2 n-1 c *n(n+l) 

The Taylor expansion for a vector function can be written 
• (d-V) B 

(16) f(d+r) - l f(r). 

n-o n! 

• 4 4 

Aj is a function of s - d*r, where r is the vector of the field point with 
respect to an origin at the center of the sphere and d is the displacement of 
the center of the sphere from the axis of the loop. Employing the Taylor 
exp)ansion, 


• (d*V) k fj lb * P 1 (n*c) 


(17) Aj - L 


{ 


A 


2 n-1 c" +1 n(n+l) 


VxP (n*r)r"r}. 


k-o k! 

The displacement vector d has been chosen so that it is perpendicular to the 
axis of the loop (d*n - 0). By the choice of the Legendre polynomial the axis 
of the loop is designated along the z-axis. It is natural then to designate d 
as being along the x axis so that 
9 


M.vr 


(d- 


9x 


■) k . 


In this case, a useful identity is [8, p. 361] 


3 ( * 1 
(18) P*r"e iBf - — 

ax 2 

Also true is, 

3 . 


(itH-n)(nH-n-l)P^ ,_ |e 1 ( "~ l ’ - P",e 

n — 1 n — 1 


+!!(■+!>♦ 


n-1 


r - e , which implies that d — r - d, so 

3x * 3x 


(19) (d*7)P^r n e iB+ r - - r n_1 r 
" 2 


(m+n) (rtH-n-l)p" *e 

n — 1 


-l^i ( m-1 ) ♦ _pm-H gi ( m+1 ) ♦ 

n-1 


+P B e iB+ r n d 

n 
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Hie important aspect here is that even after differentiation the basic form of 
each term is preserved. Each term is either of the form 
pJe lk Vr or pJe lk Vd. 

Since (d*7)d - 0, repeated application of d*V yields only more terms of 
similar form. This makes it possible to find a closed form for any power of 
d»7 operating on P n r n r. The derivation is given in Appendix C. The result 


is 

(20) (d*V) k P r n r 


d k k 


I ( k ) 


nl 


2 n^o " (n— 2m) I 


(-1) k ““p k-2 "e 1<k_2 " ) + r n-k r 


Kd 


k-1 


sk-1 


Y ( ‘- i ) JY , .!)■>— 

m-0 “ (n-2m) 1 n_k+1 


Using this result in equation (17) and the fact that d*7(7xA) ■ 7x(d*7A) we 
obtain 

n! 

kl ir^O (n-2m)l 


(21 )Aj - 7x 


' A 

p £b ® P 1 (n*c) n 


2 n-1 c n+1 n(n+l) k«0 


d 1 kn 

Z ( ) k Z (-l) k "*- 




for j r | >c . 


pk 2 m 

n - k 


(n-r)e 


i ( k-2m ) f 


k (r+d) . 


III. Calculating the total VCCtOl AUwAUC 

The differential equations for the vector potential are linear so 
solutions can be superimposed. Therefore linearly independent terms in the 
harmonic expansion for the inducing vector potential can be separately matched 
to boundary conditions to determine corresponding terms for the vector 
potential inside the sphere. The terms can be recombined to give the total 
vector potential inside the sphere. 
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As has been mentioned, all of the terms of A 1 are either of the form 
Vx[P^(n*r)e 1 "*r"r] - A h> r n rxV{p£(n-r > )e 1 " f ] 
or 

A» yx i rBp !!< n, *> eiMd ]* 

The first are recognized as vector spherical harmonics. The general solution 
to Laplace's equation (using only terms defined for all n*fc) is 

(22) A e - E E A h> [r n +B n> r -n ” 1 ][rxV(P£(n*r)e i "* )]. 

n-o m— n 

The potential outside the sphere is a combination of that due to the current 
loop and that due to the eddy currents in the sphere. The potential due to 
the eddy currents must vanish as r approaches infinity so A ni> must be 
completely determined by the current loop. There are no r -n-1 terms in the 
vector potential of the current loop so B is completely determined by the 
eddy currents. Immediately A nB is known. 

Likewise, the general solution to Helmholtz's equation can be writtin 

n-o m— n 

where 

I n (kr) - i” n J n (ikr) 

K n (kr) - l/2ni n+1 (J n (ikr)+iY n (ikr)] 

J (ikr) is a Bessel function of the first kind 

n 

J (ikr)cos(ikrn) - J (ikr) 

Y n ( ikr ) - - 

sin(ikrn) 

The function K n+1/2 is singular at the origin, which is within the sphere, so 
D -0 for all n,m. So 

n in 

(23) \ - t r - 1/2 C„.I n . (kr)rx7(^(n-t)e‘*'). 

n=o 
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The boundary conditions are met if (Appendix D) 


(24) C 


(2n+l)a 


n-l/2. 




A Einilar process is used for the terns of the form tx(pV*" + r^e^ ). 
For these terms, the curl is taken first and is broken up into its cartesian 
components. From Appendix E 

A nB Vx(p£( n ‘r) el "* rn e x ) - A nB {(iiH-n)pJ|_ 1 (n*r)e iM r n_1 e y - 


- [ (m+n)(m+n-l)P*~* (n^rje 1 <n ' 
2 


“♦ r n-1 +p" + * (n^rje 1 *' 

n — l 


l+ 1 , V -1 Je E ). 


In cartesian coordinates V 2 AeV 2 A x e x +V 2 A y e y +V 2 A x e t . Each component of 
Helmholtz' s vector equation is an identical Helmholtz scalar equation. Prom 
[p. 375) and the arguments stated for the vector spherical harmonics, the 
general solutions for A ± and A o are 


I \l(t*st.r— 1 IPjIn- V 

n-o 


j*x,y,z 


K - 1 C n l 

n«o 

The similarities between cnese equaulwiio (77) » r ** 

obvious and thus 


(2n+l)a n_1/2 A j 

n in 

ki ..i/=< k » 


The complete expression for the vector potential inside the sphere is 
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A/ 0 Ib 

(25) A. - 

* 


£A £ 

n_l" Ir.r 


Z(-l) k_ " ( k ) 


ni 


t~ 1/2 [I 


(kt ) 




. k - 2 ■ /v i ( k-2a ) f 

< VP n-k' n * r > e > Xr+<iI »-k-l/2< kr > C »-k- 1 


(n— 2m)P^:^" x (n»r)e* <k ~ 2 * , *« 


- - ((^)(n-2 n Hl)p;: k 2 ^W^)e i(k - a -- 1, %p;: k a :;U-r)e i,k ‘ 2 - +1> *e 8 ) 
2 • 


} 


where 


P 2 (n*c) (2n+l)a n-1/2 

^ c n+1 n(n+l) ' C ° kl n _ 1/2 (ka) 

IV. Calculating the average power dissipated An the sphere 

For a sinousoidally time dependent current density, the average power 
dissipated is 


P-= — J* • JdV. [3 p. 369] 

2a J 

a is the conductivity. 

The factor of one half arises from the time averaging over one period. 

tie inn emiaHnn (9) 

J * -icwA, 

the average power becomes 

(o 2 a r 

P - I A* -AdV. 

The displacement d of the center of the sphere from the axis of the loop 
is small, so only terms to first order in d will be considered. This can 
be justified by examining equation (25). Each term of A ± is propotional to 

.k n-k-1/2 , n-k-1/2 

a a da 
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The loop must be outside the sphere so a<c and d is assumed small so (4/c)is 
small. From equation (25), to first order in d 

(27) V' 0 “ "Y" n f x V" l/a lWa< kr > C B VP B ( n *r)xr 

% 

~ I n-l/2 <kr)C n-X dt7(P n-l (n '^ )COS+)Xr " ^n-l (n *^ ,e T “ P n-1 ( n *^)* in 4®. U 

Using -n(n-l)P" B _ 1 - P B-l (7, p. 560) 

e i + +e” M e 1 * -e 1 * 

and cos* ; sin| » 

2 2i 

Study of equation (27) reveals that the angular part of the integral 
(equation (26)) to first order in d has terms of the form 
(28a) | (7P n (n*r) x r) • (TP^n^r) x r)d2 


(28b) | (7P n (n*r) x r) • (y[pj_ 1 (n*r)cos4>] x r)dS 

(28c) <£ (7P B (n‘r) x r) • (1P 1 _ 1 (n*r)e y + pJ_ A (n*r)sin$ e e )dQ 

These are evaluated in Appendix F with the results 
| (7P n xr) • (y[pJ_ 1 cos*l x r)dC - 0 

£ (7P n xr) • (1 P 1 _ 1 ^ r + pJ_j sin* e E )dQ - 0, and 


£ (7P n xr) • (yPjXrJdfi - 


4nn(n+l) 

2n+l 


°n,l • 


Hius, all of the first order terms vanish and only the zeroth order 
contributes to the power dissipation. The expression for the power is now 


coo 2 u I 2 b 2 • 4nn(n+l) 

° ° j; A 2 

8 n=l " 2n+l 




Following the logic of Sonnabend (see Appendix I) 


I n+1/2 (x) — - ; r~a. 

n+1/2 (2nx) 1/2 
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80 


l Cn | 2 | I „ + l/2< kr) | 2rdr ~ 


i 3 a 2n (2n+l) 2 


where 6 ■ ( ) 1/2 is the skin depth . 


Thus, 


/i o wo 


. ow 2 // nS 3 I b 2 • 

P . 2 2 i a 2 (2n+l)n(n+l)a 2n 

4 n-1 


nl 2 b 2 - |P 1 (n-c)| 2 (2n+l)a 2n 

- E . 

oS n-1 n(n+l)C 2n+2 

V. Calculating the force on the sphere 

The force on the sphere due to the external magnetic field is 

F - J JxBdV - Re {-icuo |a £ x (VxA. )dV} 

where A t is the total vector potential inside the sphere and Aj is the vector 
potential due to the current loop alone. The interaction between the eddy 
currents and the field of the eddy currents is of no interest. The sphere 
cannot move itself and thus that force must be negated by internal stresses 
ahout which little is known. 

Evaluating tfxAj to first order in d, there are terms of the form 
(Appendix H): 
zeroth order 

(29a) 9x(Vr n P n xr) - n(n+l)r n_1 P n r - (n+l)r n_1 sin0P n ’ 0 

, ^ dP n (n*r) 

where P_ (n*r) - — 

d(n-r) 

first order 
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zeroth order 

(30a) r" 1/a C ll I B+1/I (kr)OP n xr - r- 1/J C„I„ <1/I <kr) S inep; + 
first order 

(30b) -r" 1/2 C n _ 1 I n l jtkrJVtp'cos+lxr - 

1/2 ~ S111 * 1 ~ 

- r v ♦ - — fi-.w 

sine 

(300 t- 1/! C 1 ,_ 1 I 1> . 1/I (kt)n P n .,e v 

(30d) r -1/J C n _ 1 I |l _ 1/J (kr)E^_,6in4e 8 

From the cross product of the zeroth order terms of A t with the zeroth 
order terms of VxA x the zeroth order terms of the force are obtained. From 
results in Appendix K the force to zeroth order, averaged over one period of 

*>• « *- f i e 


-//^o'b 2 !! Z ?„ (n*c)P * +1 (n*c)a 2 

n-l — 

2 2 2n+3 (n+1) 


This is identical to the force, found by Sonnabend Ip. 54], on a sphere on the 


axis of the loop. 

The first order term for the force is found by taking the cross product 
between the zeroth order terms of A t with the first order terms of A ± , and the 
cross product between the first order terms of and the zeroth order terms 
of Aj . The calculations are in Appendix J. The results is 
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-IV/ I 2 b 2 d • |P 2 (n*c)| 2 a 2n_1 (n 2 +2n-l) 

r 1 - — £ { — 

4 n-1 c 2n+2 n(n+l) 

P 2 (n*c)P 2 +2 (n*c)a 2n+1 (n+1) 
c 2n+4 (n+3) J 

Hie zeroth order force is one strictly along the axis of the loop, and the 
first order force is one strictly in the direction of the displacement from 
the loop. At first glance the first order force appears to return the sphere 
to the axis of the loop, however evaluation of the summation should be done to 
be sure. The first term in the summation is definitely restoring because 

|p 2 (n*c) | 2 >0 

2 n - 1 

8 

and it should dominate because of its dependence versus the other terms 

p2n + l r 

a 2n+2 U 

a 

dependence. 

Jn+4 

C 

VI Calculating the torque about the center of the sphere 
The torque is found from 

r - jr X (J X B)dV. 


J x B was calculated to first order in d when calculating the force. Finding 
the terms of r x (J x B) is then trivial. Calculating the torque then becomes 
simply a matter of trudging through the integrals as was done in calculating 
the force. This is done in Appendix L. The results are that there is no 
torque to zero order, and to first order the torque is 


r 1 _ 

IV/ I 2 b 2 d 
0 

• fp 2 (n*c)a 2n-2 

P 2 (n*c)P 2 + 1 (n-c)a 2n+1 ( 2n 3 +3n 2 +3n-2 ) ' 


"** 4 

n-ll c 2n+1 

c 2n+3 (n+l)(n+2)(2n+l) 


Thus the axis of rotation is perpendicular to the displacement off axis. For 
a displacement in the opposite direction (d ■+ -d) the torque changes sign. 
This means that as long as the sphere stays near the axis of the loop, the 
angular acceleration will average to zero. 
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APPTODIX A 


From equation (13) the vector potential of the loop is 


Expand 


1 


1*1 



in Legendre functions. 


* 


1 • s 

* I P Cb-c), ' I7,p. 539] 

|r| n-o c n+ 

c is the vector to the loop, s to the field point. Functions of c can be 
separated from functions of s by using the addition theorem for Legendre 
functions in terms of spherical harmonics. Let n be the normal unit vector to 
the loop define the z-axis. 


P n (s*c) 


~ ^ n (n-m)! 

P (n*s)P (n*c)+2 Z 

n " m=l (n+m)l 


p£(n*s)p£(n*c)cosm(f c -^ g ) 


where P* (x) 

n 


(l-x 2 ) 


1/2 


d" 

dx" 


P n (x). 


[7, p.582] . 


For the purpose of simplicty, the coordinate system can be defined where ^* 0 . 
Since c«|c| is constant over the loop as in n*c, the line integral amounts to 
an integration over <p c . 


n2TT 


d> dl 


bd4>(-sin<f<e +cos<f>e ) 

x y 


Using the following orthogonality relations, 


sin* cos in$ d$ 

OTT 

cos$ cosm<|> d<|> « 

Jo 



* o 


JIS. 


m , 1 



n^l 


[7, p. 435] 
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p lb «* P* (n*c) 

A, - I — sV(n.£)e. 

2 n-1 c n+1 n(m-l) T 

Hie coordinates were fixed relative to s when was set to zero. This 
can be "unfixfed" by letting e y - \ , and the coordinates are again,, independent 
of 8. Thus 


p o Ib 


n-1 c" +1 n(n+l) 


P* (n*c)P* (n*s)f s . 
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APPENDIX B 


TO see that s n P 1 (n*s)* - 7 fP(n > s)s n s] ; first note that 

n Bin 

0 

7(s n P (n*s) ) - ns n_1 P„ (n*s)s + s n_l P(n*s)0 . 

* 89 , 

From the chain rule, 


36 . 


P„(n-s) 


3 

3(n-s) 


* A . 

P n (n-s) 


3(n*s) 

30, 


- -p ’ (n*s)sin0 - -P*(n*s). 

n Bn 


sxs-o and sx0,«s+, , so 
7 (s n P (n*s)) x s > 

s n 

7,x{P n (n-s)s n s] 
and 7, xs - 0. So, 

V x [P (n*s)s n sj ■ 

b n 


* (”•*)♦, 

■ V [P (n*s)s B ]xs + P (n’S)s B P xs 

Bn n S 

7, [P n (n*s)s n ]xs « s n P*(n*s)+,. 
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APPENDIX C 


Using equation (19) with in-o, 


(d*V) (P r n r)» — r"' 1 (n(n-l)P 'Je'—P,, Je 1+ ]r + P r n d. 

n ^ 2 n “ a n~i n 


Repeating this yields 


(d*V) 2 (P B *r n r) - — r n - 2 [n(n-l)(n-2)(n-3)P n :2e" 2if -2n(n-l)P*_ 2 + P n _ 2 e 2lt ]r 
4 




+ r”" 1 (n(n-l)P^ 1 e" 14 -P B ^ 1 e i4 ]d 


(— )V ' 2 l C) 

2 ir=o ( n— 2m) l 


(-l) 2 '*P 2_2 *e i<2 ' 2 " >+ 


d 1 , nt 

+ 2 ( ) I ( 2 ) 

2 m-o “ (n-2m) ! 


, 2-2»+l n-2+l. 

(- 1 ) P n - 2 + 1 r d * 


The pattern holds for further differentiations with the result 


(d-y) k (P n r n r) 


d k . n! 

— z ( ) 

4 m=o * (n-2m)l 


< k-2» ) ♦ ^ji-k f 

n — k 


kd^ 1 k x 
E ( ) 


0 k-l 

2 m=o 


n! 


(n-2m) I 


. . vk-m-l k-2»-l n-k + 1 . 

*n-k + l C “* 
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APPENDIX D 


Taking each term independently, \ | r-> - A o | r _ a implies that 

(Dl) [\.(t”+B„.r-”- 1 )rxro; e 1 -‘| r ...|C„.r 1/i I 1> . l/2 (kc)n t <e I **| r ... 

rx7P"e i -" + can be cancelled so, 

n • 

(D2) -c„a-‘ /J i > , 1/1 <l«i). 

The next boundary condition is that 7xA i | r _ a » VxA o |r«a. Taking the curl of 
equation (Dl), cancelling the angular dependence and using equation (F2) 
yields 

Carrying out this differentiation in spherical coordinates and equating 
components yields. 

(D3) * h .(na"-(n + l)B„.a-”-‘) - C n .{-l/2a !/! I n . 1/J (ka) + a- 1/2 kI- 1 , 41/! (ka)). 

From (3, p. 98] 

n+1/2 

« w > I '.../ ! "'a) - V 1/2 (ka) - — I„, 1/: (ka). 


Substituting this into (D3), solving for B nB in (D2) and then solving (D3) for 
C in terms of A yields "" 

n in n m ■* 

A a n-1/2 (2n+l) 
kl„_ w ,(ka) * 
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Appendix E 


Vx(r D £^e 1,f e a ) (r n pV" 4 )e (r'Ve 1 "* )e s 

3z 3y 

From {8, p. 361] 

r B P p e i,f - — I(BM-n)(iiH-n-l)p""Jr , '" 1 e 1( "“ 1, *+P* + Jr n_1 e i<1,+1)f ]. 

3y “ 2 

* * 

J_ r l *P^e i *^coser ,, -Ve 1 * 4 +sin 2 eP^ , r n - 1 e i * 4 

n n n 


- (mxn)!^^ r n-1 e i ” + . (using equation (L17)) 

Thus, 


«-l n - 1 i (»-!)♦ 


Vx(r n pV" + e a )-(nH-n)p£_ 1 r n_1 e i "*e [ (nH-n)(nH-n-l)P^J r n_1 e 


+p"**r n_l e i lm+i) * ]e a . 
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Appendix F 


To evaluate (28a) and (28b), refer to the theory of vector spherical 

harmonies which are defined in [9, p. 211-214] as 

-i (2n+l)(n-m) 1 

-7 " ™7{(-l)"[ ) 1/2 l*e im *} . 


[n(n+l)] 1/2 


4n(n+m) l 


Thus, 


VP n xr - -i [- 


4itn(n+l) 


] 1/2 ^ 


2n+l 


and 


yp x _Jcos4 xr - l/2V(P 1 _Je i ii) iii) *-l(l+l)pJ_ l e" if ] xr 
-i 4 it 

( ) 1/2 1(1-1)(X. .-Xi .). 

2 21-1 ^' l ' 

The vector spherical Harmonies satisfy the orthogonality relation 


% 


[9, p. 211] 


It can be immediately seen that (28b) vanishes, to evaluate (28a), 


[ (TO^rMVPj 


xr)d2 - — 4n[- 


n(n+l) 


2n+l 


l 1/2 l- 


1 ( 1 + 1 ) 


21+1 


•J 1/2 fxno-X 1( 


dQ 


4n[ - 


n(n+l) 

2n+l 


] 1/2 [ 


1 ( 1 + 1 ) 

21+1 




X *x dfi 

A no A lo u 


4nn(n+l) 

2n+l 


n 1 


TO evaluate (28c), the following relations are useful: 


9P 


i) VP xr «= - 

n 


99 


ii) ♦ • e y = cos4> 

iii) ♦ • e t - 0 
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cos$d4> - 0 



Using i), ii) and iii), (28c) becomes 

| (VP'XrM^lF^ + • B Pj_ 1 sinfldfi - j - lPj.^os^da 

ft 

which, using iv) vanishes. 
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APPENDIX G 


The following argument is from [1, p. 53]. 

(e*-(-l) n e -x (n+1) ! 




(2nx)" 1/2 I 


l-o lKn-UM-fe) 1 

It has already been assumed (assumption 5 of the introduction) that the skin 
depth & is small compared to the radius of the sphere. This means that the 

* 4 

current density will be negligible except when r~a. In the equation for the 
current density the argument 

r r 

x-kr-(i// «o) 1/2 r - (2i) 1/2 — - (1+i) — . 

8 8 

But for r~a x ~ (1+i) a/8 is a large number, so the summation e x dominates 
over e“ x . (-2x) -1 goes quickly to zero as 1 increases, so the l-o term of the 
summation dominates . Thus 

\ 

Remembering that most of the current density is at r~a and that over small 
variations in r, I_, ,,(kr) will fluctuate much more than small powers of r, it 
is reasonable to let all of the functional variation lie in the exponential 
and set r-a in the radial integral for small powers of r. Thus 

Un+ij a 


n+l/2 


(kr) ~ (2nkr) -1/2 e k r 


fa 

l C n 

Jo 


I IWa<M 


fa 

| 2 rdr - 

Jo 


k 2 e 2ka 


e 2kr dr 


O G \ < a / 
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APPENDIX H 


Zeroth order terms of \ have the form 7(r n P n )xr, and first order terms are of 

the form 7(r n P )x e and 7(r n_ 1 P* , cos*)xr. So the terms of 7xA. have the 
form 

7x l7(r* l P B )xr], Vx[7(r n- 1 P n ^ 1 cos*)xr] 

and 

7x[7(r n P n )x e E J. 

For a function f(r) Y (©,♦) where V 2 [f(r)Y(0,4)]-o 

7x[7(fY)xr] - <r-7)7(fY) - [7(fY)*7]r + (7*r)7(fY) - [V*(ft)]r 
- (r-7)7(fY) + 27(fY). 

For f (r) - r n , 

n 1 - n-X 3Y - r "" 1 3Y - 

7(fY) - nr n Yr + r — © + — ♦ 

3© sin 0 3* 

3 Y r n_1 3 y 

and ( r*7)7(fY) - n(n-l)r n_1 Yr + (n-l)r "' 1 — 0 + (n-1) — f. 

30 sin© 3<t> 

Using these results, 

Vx[7(r P n )xrj - n(n+l)r n- 1 P n r - {n+l)r n_1 sin© P n » 0 

and 

7xf 7( r n_ 1 P* ,cos<|>)xr] = n(n-l)r" - 2 P B _ 1 cos<|>r-nr n_ 2 sin0P B j 1 cos«|>0 
r n_ 2 sin+ ^ 
sin© 

similarly, 

7x(7(fY)x e^ ] - (e x -7)7(fY) - l7(fY)*7]e s + [7*ej7fY - 7 2 (fY)e i 
3 

7(fY). 

3x 

so, 

3 , - sin<|> ^ 

7x[7(r n P )xe ] - — 7(r n P) - r n ' 2 [-(n-l)P cos+r + sin© P' © + P ' ♦). 

n * 3 x sin© 
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APPENDIX I 


Using (29a) and (30a), the zeroth order terms of ^ x (VxAj) have the form, 
r“ 1/2 C I .. (kr)r 1-1 (1(1+1 )P. sinQP,, f #xr-( 1+1 )sineP/sin6P 'fx©] . 

n n+i/^ i n * n 

4xr - $ - -sin0e s + cos© cos*e, + cos© sinfe . 

% 

-fx© - r - cos ©e B + sin© cos+e, + sin© sin|e y . 

Performing the angular part of the integral. 

1 

Use sin© p ' - P*, sin© P x (P l+ | - P X _J J (from Appendix N) 

21+1 


I. 


2 * 




and that sin$d$ 
o 


cos$d4> »o. Then 


:il j 


(II) l(l+l)P 1 sin©P n '©d2 


-4nn(n+l) 

2n+l 


'n(n-l) 


2n-l 


n , 1 ♦ 1 


(n+1) (n+2) 
2n+3 


n ,1-1 


Then use sin©P 1 • « pj , sin ©P^ - P* , and 


cos©p| - [(1+1)P 1 \ + 1pJ + 1 ] to find 

21+1 


(12) |sin©P x r sin©P n r l(l+l)rdS 


4nn(n+l) 

2n+l 


(n+2)' 


2n+3 


, 1 - 1 + 


n(n-l) 


2n-l 


n , 1 + 1 


Adding un ana \i.t) yielJ* vV '“ 


Jl(l+l)P 1 sin©P n '©+ (1+1) sin©P 1 ' sin© P n ' r 
Using the results from Appendix G, the radial part is 


4itn(n+l) (n+2) 


2n+l 


e « , 1 + 1 




5(2n+l)a 


2n+ 1 


(l+i)k 

Putting in the necessary constants and averaging over one cycle yields, 

-fj I 2 b 2 n 00 P*(n*c)P* (n*c)a 2n+1 

a o o _ n n + i 


i 

n=l 


c 2n+3 (n+1) 
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APPENDIX J 


The cross product between (29a) and (30b) yields 


(Jl) r" I/2 C 0 _ l I n _ 1/2 (kr)r 1-1 { 



r-sine e B I 

a) -sinep 2 '.cos+Kl+DP, 

|cos6 cost e i I 

n - jl x 

[cos© 6int e y J 

; sint 

r-sint el 

b) P^_ 1(1+1 )P X 


sine n 1 

[cost e y J 


("cose e s 1 

c ) -sin 2 6P 2 ' . cost( 1+1 )P, ' 

sine cost e K 

n — i x 

[sine sint e J 


Similarly for (29a) with (30c): 
<J2> c” 1/2 C n _ 1 I n _ 1/3 (kr)r 1-1 ( 

a) nP n _ 1 l(l+l)P 1 


cos© e x 
-sin© cost e B 


b ) -nF n _ t < 1+1 ) sir.6P 1 ' 


-sine e K 
-cose cost e t 


For (29a) with (30d) 

<J3) I 


a) p^jSintKl+DPj 

b ) -P 2 _ x sin6( 1+1 ) sinOP' 
For (29b) and (29c) with (30a) 


i 

7 

e.. 


sine cost e 
-sine sint 
Tcose cost e, 
-cose sint 




(J4) t- 1/2 C n I I ,, 1/3 (kt)t 1 - ! 


a)(l-l) 2 P 1 l; 1 cost sin6P n ' 


-sine e s 
cose cost e K 
.cose sint e 
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b)sin J 0p| ' (l-l)P n 'cos4> 


cos 9 e s 
sin0 cos4> e x 
cos© sinf e J 


From the ♦ dependence, it can readily be seen that only the x-components of 
these terms will survive integration. 

% 

Taking (Jl) first and completing the angular integration, use the 
following identities: 

i) Jcos 2 fd+ - Jsin 2 f34> - n 


ii) -sinGcosGP 1 ’ , + P 1 , * n(n-l)P„ -(n-l)sin6P 2 . 

n — l , n — l n — z n — j. 

sin@ 

iii) sin 2 0P 2 .-i~ [-(n-l) 2 P 2 +n 2 P 2 _ 2 ] 

2n-l 

iv) sin0P 1 * - pj . 

Then, 

sin 2 $ 

f [ -sinGP 1 ',1(1+1) cos© cos 2 4>P, + P 1 ,1(1+1) P, -sin 3 6P 2 ' . ( 1+1 ) P. ' cos 2 ♦ ] da 

n- l l ,n—i x n- x x 

J sine 


.Jif'l ( n-1 )nP n _ 2 1 ( 1+1 ) - ( n-1 ) sin6P 2 _ x 1 ( 1+1 ) P x [ - ( n-1 ) 2 P 2 +n 2 P 2 _ 2 ] ( 1+1 ) pJ Idcose 

i\ 2n-l 


Use sine P 2 (P 1 + 2 - P 2 2 2 ] and this equals 


21+1 


1 ( 1 + 1 ) 


*|'n(n-l)P, | _ 2 1(1+1 )P, - (n-1)?;., lP,ii - Pili 1 

-i ZItI 


2n-l 


l-(n-l) 2 n-l) 2 P 2 + n 2 P n 2 2 )(l+l)pJ dcose. 


Using the orthogonality relation for Legendre polynomials, this simply equals 


2nn(n+l) (n-l)‘ 
2n-l 


D - • 

n , 1 
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TO do the same for (J2) use, 

1 

i) cose P x [lPj x + (1+1)P 1 + 1 ] 

21+1 

1 

ii) sin 2 0P, ' [-1(1+1)P +1(1+1)P ] 

21+1 

and orthogonality to find that 

j(nP B _ 1 l(i+l)P 1 cos0 + nP n _ x (l+Dsi^ePj ' )da 

4nn 2 (n+l)5 , 

2n-l 

For (J3) use 

i) sine P, - — [Pji, - p.lj 
21+1 

1 

ii) sine cosGPj ' - t(l+l)P * + IP 2 ] 

21+1 

to find that 

j ( -P* _ 1 sinei ( 1+1 ) r 1 sin 2 $ + P n ^ sin© cos0( 1+1 ) P x ' sin 2 ♦ )dQ 
4nn(n-l) (n+l)8 . 

n , i 

■ . 

^il~ A 

for (J4) use 

i) sine cose P ' - ((n+l)P 2 + nP 2 J 

2n+l 

ii) sinOP ' - P 1 

n n 

iii) sin 2 epj ' — [-(1-1) 2 P 2 + 1 2 pJ 2 ] 

21-1 
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to find that 

[((l-l) 2 pJ_ 1 cos0sin0P i [cos 2 + + sin s +pj' (l-l)P n ' cos 2 4»]d2 

J 

2nn(n+l) 2 (n+2) 


2n+l 


“n , 1 - 2 * 


(jl), (J2), ( J3) all have the same radial terms and constants in the total 
force expression, thus their angular terms can be added together. . 


fn(n+l)(n-l) 2 2n 2 (n+l) 2n(n-l)(n+l) 

2it + + 


2n-l 


2n-l 


2n-l 


n , 1 


2n(n + 2n-l) n (n+1) 
2n-l 


°n,l * 


Because of the Kronecker delta, the radial part has the form 


r 

Jo 


2n-l 


a 2 "" 1 exp 


f(l+i)(r-a) 


dr 


8 2 (2n-l)a 2n_1 


2i 


The radial part of (J4) becomes 


‘a 2n+l 


- k 


a 2n+1 exp 


( 1+i ) ( r-a ) 


dr 


S 2 (2n+l)a 2n+1 
2i 


Putting all of this together with the appropriate constants and averaging over 
one cycle yields the first order term of the force. 

r|F 1 /(n-c)|V'" 1 (n 1 +2i>-l) <n*c>P n (n*8) (n+1) ] 


-KP.I.'b'd 


IVI 


£ 

n-1 


c 2n+2 n(n+l) 


c 2n-f4 (n+3) 


e «- 
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APPENDIX K 


For the zeroth order part of the torque, take the cross product of r with 
(Jl) and (J2). For (J2), rxr - o, so that integral vanishes. For (Jl), rx© - 4> 
so, • 

(Kl) r x jld+DPjSineP^eda - r j 1 ( 1+1 ) P A sin©P* fdfi . 

But, $ - -sin+ e s + cos+ e y 
and J ! "sin4dt - £* cos$d+ - o. 

so (Kl) vanishes also. Thus there is no torque to zeroth order in d. 

The first order terms can be found by taking the cross product with r and 
equations (Jl), (J2), (J3) and (J4). Noting that 
rxr - o 

a 

rx© - r+ 

•A A 

rxf - -r0 

rxe x - -sin© sin$ e s + cos© e y 
rxe y « -cos© e x + sinG CGSy e g 
rxe s ■ -sin© cos$ e y + sin© sin<|> e x 

These terms are readily found. The cross product of r with equation (Jl) 
yields 

<K2) r" 1/I C n _ 1 I n _ l/J (kr)r 1 { 

[-sin* e 

a ) -sin©P* l x cos4»l ( 1+1 ) P A 

lcos+ e y 

sint [-sin© e 

b) P^_ 1 1(1+1)P 1 cos© cos$ e x 

sin© Lcos© sin$ e y 

The cross product with equation (J2) yields 
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<K3) r-‘ /2 C n . 1 I n . 1/! (kr)r 1 { 


a) nP B _ 1 1(1+1 )Pj 


-sin 2 © cos* sin* e 
cos 2 e + sin 2 © cos 2 * e 
L-cos© sin© sin* e 


b)-nP B _ 1 ( 1+1 ) sinQPJ 


-cos© sin© sin 2 * e 
sin 2 © sin* e s 
.cos 2 © cos* e 


For (J3) 


a) p^jSin+Kl+DPj 


-cos© sin© cos* e x 
-cos© sin© sin* e 
.sin 2 © e_ 


b) -E^_ 1 sin*(l+l)sinGP' 


-cos © cos* e B 
-cos 2 © sin* e 
Lsin© cos© e B 


Finally, for (J4) ' 

(K5) t- 1/: C„I n . 1/I (kr)r‘- , { 


(l-l) 2 P 1 _Jcos* sinQP^ 


-sin* e B 
Lcos* e„ 


}• 


Using 


j 

0 * sin*d* « cos*d* * cos*sin*d* « cos *sin*d* - sin *cos*d* 

Jo Jo Jo Jo Jo 


~“ti“ cv" f CTVk fhp anoular 


oruy tne y-eoniponeuLt> &ulvavc - 
part of the integral for terms of the form of equation (K2) use 


p2T\ <« y 

i) | cos *d* - n - J sin *d* 


cos© 

ii) sine - — P„_i - - nln-llf.., 

sin© 


, * 

iii) [P. ,P,dS -= 8 

J n_1 1 2n-l n 1,1 


32 


Then integration of the angular part of (K2) yields 


2it 

2n-l 


n 2 (n-1) 2 S, 


n - 1 , X 


To evaluate the angular part of (K3) use 


r ‘ in!+ ■ r 


cos ♦ - n 


* 


to eliminate the ♦ dependence. The angular part being integrated is then 
2imP B _ 1 l(l+l)P 1 cos 2 e + nnP n _ l l(l+l)P l sin 2 e + (l+DP/ Sin 2 e cose. 

Using 

cosGP' ■ 1P X + P 1 _' this becomes 

2nnP n _ 1 l(l+l)P x + nnP B _ 1 (l+l)P 1 ^ 1 sin 2 e. 

The first term can immediately be evaluted. It's integral is 

, 4im 2 (n-l) 

2wtf> , 1(1+1 )P, dcose - S . . 

h — 1 * A a n - i I 1 


Using 

, Kl-l) 

sin 2 ©P' [P._ 2 -P.] # 

21-1 


the second term yields 

2nn 2 (n+l)(n+2) 2nn 2 (n-l)(n-2) 

( 2n+l ) ( 2n-l ) n ' 1 " 1 (2n-3)(2n-l) B,1+l 

Doing the same with (K4). Use 

i) P^ - -lPj+cosePJ 

ii) cos0P 2 - [ (n-l)P n 2 +nP B 2 2 ] 

2n-l 

iii) sine P^ - P X 2 X . 

Integrating, the angular part yields 

2nn(n+2)(n-l)(n+l) 2nn 2 (n-l)(n-2) 

( 2n-l ) ( 2n+l ) "' 1_1 (2n-l)(2n-3) n,1+1 
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Equations (K2), (K3) , (K4) all have the same radial dependence so their 
angular parts can be added. This yields 


2nn(n+l)(n+2) 


2n+l 


^n,X-X + 


2im (n-l)(n+l) 
2n-l 


«,x*x * 


For (M5) use 


sinW • m p . 

n - n 


The integral then yields 


2nn (n+1) 
2n+l 


°n , 1 - X * 


Integrating the radial part of these terms is similar to that integrations 
done in finding the force. 

f C »-l I »-l/ 2 < kr)rl+2 " 1/2dr - J C »-l I n-l /2 < kr > r " +3 " 1/2 dr f ° r *n,X-l terttS 
“ f C „-l I n-l/ 2 < kr > r ' , + 1 ~ 1/2 dr fOC *».X*1 tennS ‘ 


and J C n I n+1/2 (kr)r 1+1 1/2 dr - J C n I n+1/2 (kr)r n+2 1/2 dr for (M5) terms. 

Again, using the same approximations as were used in finding the force and 
power, 

\» 2n+1 


(lti)k 




1/2 


S(2n-l)a 


2 n-X 


(l+i)k 




S(2n+l)a 


2 n + 1 


(l+i)k 

So, the final expression for the first order term of the torque is 


- e y^o I o 2b2d Z 

r - — z 

4 n=l 


P* (n*c)P^_ 1 (n*c)a 2n_1 P* (n-c)P* +1 (n-c)a 2n+1 ( 2n 3 +3n 2 +3n-2 ) 
+ 


2 n + 1 


c 2n+3 (n+1) (n+2)2n+l) 
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APPENDIX L 


From [7, p. 541] 

(LI) P B+ J(x) • (n+l)P n (x) + xP n '(x) 

(L2) p a_i<x) - -nP B (x) + xP n '(x) * 

(L3) (l-x 2 )P n '(x) - nP n _j(x) - nxP n (x) 

(L4) (l-x a )? n '(x) - (n+l)xP n (x) - (n+l)P B+1 (x) 

From (3,p. 162-3) 

L5) P ’II< X > * (n+®+l)(l-x 2 ) 1/2 p£(x) + xP^ +1 (x) 

<L6) p"J<x) ~ C!< x > * (2n+l)(l-x 2 ) 1/2 P^(x) 

(L7) p"J(x) - (m-n)(l-x 2 ) 1/2 P^(x) + xP^ +l (x) 

(L8) P * +1 (x) - 2mx(l-x 2 )~ 1/ V(x) + (nw-n) (n-nH-l)P^ -1 (x) - 0 
(L9) (m-n-l)pjj +1 (x) + (2n+l)xP^(x) - (m+n)!^ (x) - 0 
(L10) (l-x 2 ) 1/2 P^' (x) - -mx(l-x 2 f 1/ V(x) + P£ +1 (x) 

(Lll) - -l/2(nH-n)(n-mfl)P^ _1 (x) - 1/2 p£ +1 (x) 

(L12) - mx(l-x 2 )“ 1/2 P^(x) - (n+m)(n-nH-l)P^ _i (x) 

(L13) 2m(l-x 2 )“ 1/2 p£(x) - P*!j(x) + (m+n-l)(nH-n)p"~* (x) 

(L14) (l-x 2 )P^'(x) - mxP™ (x) - (irw-n) (n-m+l) (2n+l) -1 (p£ + 1 (x) - Fj_ t (x)] 

(L15) - (n+1) x Pj(x) - (n-nH-l)p£ +1 (x) 

(L16) - (2n+l) -1 [ (m-n-l)np£ +1 (x) + (n+1) (m+n)P^_ 1 (x) ] 

(L17) - -nxpjj(x) + (nH-n)I^_ 1 (x) 
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APPENDIX M 


Another expression for the current density in the sphere 

Ludwig Hannakam found an exact integral form for the eddy current density 

in conducting .solid sphere in the presence of an alternating current loop of 

% 

arbitrary shape or orientation to the sphere. lie first constructed a surface 

by drawing lines through the center of the sphere and the current loop. He 

then showed that the vector potential due to the loop can be calculated by 

considering radial magnetic dipoles on the constructed surface. The problem 

of a conducting sphere in the presence of a magnetic dipole of oscillating 

dipole moment has a known solution. He partially completes the integral over 

the surface leaving only an integral over the current loop. His solution is 

iu#/ 0 l • 2n+l 1 _ 

J - I f (kr)P(r) 

4n n-1 n(n+l) a n 


where a - ( l)n f (ka) + kaf , (ka) 

n n n-i 


o 

2 


f„<kr) - (— ) ,/J I.. 1/1 (kt) 
kr 


a 

— 'V r . ,n , 

r - I. * y V , 

J r 

c 


dS_ 


c v- .,0.2 X . 0 r ,0 

"n ' c ' c * c 


dS. 


C ^ A A 

— WrlP "fr»r )1 


a is the radius of the sphere 

r e is a vectbr to a point on the current loop (source point) 

r is the field point in the sphere 

dS c is a line element on the current loop 


I ■ I e iwt is the current in the loop 


36 


REFI2U5NCES 

1. Sonnabend, D., "Study of a New Drag-Free Air Density Instrument Requiring 
no Propellant" Stanford Univerrsity 4-1976. 

2. Sonnabend, D. , "SemiDrag Free Gravity Gradiometry”, J. Astr. Sci. 33, 

353-365 (1985). — 

3. Smythe, W.R., Static and Dynamic Electricity, McGraw-Hill, New York 
(1968). 

4. Tegopoulo? , J.A. and Kriezis, E.E., Eddy Currents in Linear Conducting 
Media , Elsevien, N.Y. (1985). 

5. Hannakan L. "Wirbelstrome in der Kiegel bei Beliebig Geformter Er regender 
Leiterschleife." Z.F. Angewandte Physik 32, Heft 5-6 (1972). 

6. Miller (Jr.) w. , Symmetry and Separation of Variables, Addison-wesley, 
Reading (1977). 

7. Arfken, G., Mathematical Methods for Physicists, Academic Press, N.Y. 
(1970). 

8. Bateman, H., Partial Differential Equations of Mathematical Physics , 
Cambridge, University Press (1964). 

9. Jackson, J.D., Classical Electrodynamics , Wiley and Sons, N. Y. (1962). 


9 


37 



